WHEN DOES A SCHRODINGER HEAT EQUATION PERMIT 

POSITIVE SOLUTIONS 
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Abstract. We introduce some new classes of time dependent functions whose denning 
properties take into account of oscillations around singularities. We study properties of 
solutions to the heat equation with coefficients in these classes which are much more 
singular than those allowed under the current theory. In the case of L 2 potentials and 
L 2 solutions, we give a characterization of potentials which allow the Schrodinger heat 
equation to have a positive solution. This provides a new result on the long running 
problem of identifying potentials permitting a positive solution to the Schrodinger equation. 

We also establish a nearly necessary and sufficient condition on certain sign changing 
potentials such that the corresponding heat kernel has Gaussian upper and lower bound. 

Some applications to the Navier-Stokes equations are given. In particular, we derive a 
new type of a priori estimate for solutions of Navier-Stokes equations. The point is that 
the gap between this estimate and a sufficient condition for all time smoothness of the 
solution is logarithmic . 
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1. Introduction 

In the first part of the paper we would like to study the heat equation with a singular 
L 2 potential V = V(x,t) , i.e. 



Am + Vu - u t = 0, 

u(x, 0) = Uq(x) 



in R n x(0,oo), V G L 2 (R n x (0, oo)), n>3, 
x e R n , M <= L 2 (R n ). 



Since we are only concerned with local regularity issue in this paper, we will always assume 
that V is zero outside of a cylinder in space time: -8(0, Rq) x [0, To], unless stated otherwise. 
Here Rq and To are fixed positive number. The L 2 condition on the potential V is modeled 
after the three dimensional vorticity equation derived from the Navier-Stokes equation. 
There the potential V is in fact the gradient of the velocity which is known to be in L 2 . 
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The unknown function u in (jl.lj) corresponds to the vorticity which is also known to be a 
L 2 function. 

We will use the following definition of weak solutions. 

Definition 1.1. Let T > 0. We say that u e L} oc (R n x (0, ,T)) is a solution to (HHJ) if 
V(.)u(-,t) € Ll c (R n x (0,T)) and 

f u (x)(J)(x)dx + / / u(f) t dxdt + / / uA(f)dxdt + / / Vu<f)dxdt = 

R n JO JR" JO ifl JO JR" 

for all smooth, compactly supported vanishing on R n x {T}. 

It is well known that L 2 potentials in general are too singular to allow weak solutions 
of (jl.lj) to be bounded or unique. Therefor further assumptions must be imposed in order 
to establish a regularity theory. The classical condition on the potential V for Holder 
continuity and uniqueness of weak solutions is that V G L v f c with - + ~ < 2. This 
condition is sharp in general since one can easily construct a counter example. For instance 
for V = a/\x\ 2 with a > 0, then there is no bounded positive solution to (jl.lj) (c.f. [BG]). 
In fact in that paper, it is shown that if a is sufficiently large, then even weak positive 
solutions can not exist. There is a long history of finding larger class of potentials such 
that some regularity of the weak solutions is possible. Among them is the Kato class, 
time independent or otherwise. Roughly speaking a function is in a Kato type class if the 
convolution of the absolute value of the function and the fundamental solution of Laplace 
or the heat equation is bounded. This class of functions are moderately more general than 
the standard L p,q class. However, it is still far from enough for applications in such places 
as the vorticity equation mentioned above. We refer the reader to the papers [AS], [Si], 
[Z], [LS] and reference therein for results in this direction. The main results there is the 
continuity of weak solutions with potentials in the Kato class. In addition, equation (jl.lj) 
with V in Morrey or Besov classes are also studied. However, these classes are essentially 
logarithmic improvements over the standard L p,q class. In the paper [St], K. Sturm proved 
Gaussian upper and lower bound for the fundamental solution when the potential belongs 
to a class of time independent, singular oscillating functions. His condition is on the L l 
bound of the fundamental solution of a slightly "larger" potential. 

In this paper we introduce a new class of time dependent potentials which can be written 
as a nonlinear combination of derivatives of a function. The general idea of studying elliptic 
and parabolic equations with potentials as the spatial derivative of some functions is not 
new. This has been used in the classical books [LSU], [GT] and [Lieb]. Here we also allow 
the appearance of time derivative which can not be dominated by the Laplace operator. 
Another innovation is the use of a suitable combination of derivatives. The class we are 
going to define in section 2 essentially characterize all L 2 potentials which allow (jl.lj) to 
have positive L 2 solutions. 

The question of whether the Laplace or the heat equation with a potential possesses a 
positive solution has been a long standing one. For the Laplace equation, when the potential 
has only mild singularity, i.e. in the Kato class, a satisfactory answer can be found in the 
Allegretto-Piepenbrink theory. See Theorem C.8.1 in the survey paper [Si]. Brezis and 
J. L. Lions (see [BG] pl22) asked when (jl.lj) with more singular potential has a positive 
solutions. This problem was solved in [BG] when V = a/\x\ 2 with a > 0. In the case of 
general time independent potentials V > 0, it was solved in [CM] and later generalized in 
[GZ]. However the case of time dependent or sign changing potentials is completely open. 
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One of the main results of the paper (Theorem 2.1) gives a solution of the problem with L 2 
potentials. The main advantage of the new class of potentials is that it correctly captures 
the cancelation effect of sign changing functions. Moreover, we show in Theorem 2.2-3 
below that, if we just narrow the class a little, then the weak fundamental solutions not 
only exist but also have Gaussian upper bound. A Gaussian lower bound is also established 
under further but necessary restrictions. 

Some of the results of the paper can be generalized beyond L 2 potentials. However we 
will not seek full generalization this time. 

Before proceeding further let us fix some notations and symbols, to which will refer the 
reader going over the rest of the paper. 

Notations. We will use R + to denote (0, oo). The letter C, c with or without index 
will denote generic positive constants whose value may change from line to line, unless 
specified otherwise. When we say a time dependent function is in L 2 we mean its square is 
integrable in R n x R + . We use Gy to denote the fundamental solution of if it exists. 
Please see the next section for its existence and uniqueness. The symbol Go will denote the 
fundamental solution of the heat equation free of potentials. Give b > 0, we will use <jr& to 
denote a Gaussian with b as the exponential parameter, i.e. 

^ = ^ l t;y >a ) = ^^-^e-^l*/(*--). 
Given a L] oc function / in space time, we will use gt, * f(x, t) to denote 

/ / 9b(x,t;y,s)f(y,s)dyds. 

When we say that Gy has Gaussian upper bound, we mean that exists b > and c > 
such that Gy(x, t; y, s) < cgb(x, t; y, s). The same goes for the Gaussian lower bound. 

When we say a function is a positive solution to (jl.l)) we mean it is a nonnegative weak 
solution which is not identically zero. 

Here is the plan of the paper. In the next section we provide the definitions, statements 
and proofs of the main results. In section 3, we define another class of singular potentials, 
called heat bounded class. Some applications to the Navier-Stokes equation is given in 
Section 4. 

2. SINGULAR POTENTIALS AS COMBINATIONS OF DERIVATIVES 

2.1. Definitions, Statements of Theorems. 



Definition 2.1. Given two functions V £ L 2 (R" x R + ) ; / £ £* oc (R n x R + ) and a > 0, 

we say that 

V = Af-a\Vf\ 2 -f u 

if there exists sequences of functions {Vi} and {fi} such that the following conditions hold 
for all i = 1,2...: 

(i) . V £ L 2 (R n x R+), Afi £ L 2 (R n x R+), d t fi £ L 2 (K n x R+), £ £ C(R n x R+). 

(ii) . Vi -> V strongly in L 2 (K n x R+), \Vi\ < \V i+1 \, fo^f a.e. and fi{x,0) = f i+x (x,0). 
(Hi). 

V^Afc-alVf^-dtU 
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Here we remark that we do not assume Af, |V/| 2 or dtf are in L 2 individually. This 
explains the lengthy appearance of the definition. 

The main results of Section 2 are the next three theorems. The first one states a necessary 
and sufficient condition such that (jl.lj) possesses a positive solution. The second theorem 
establishes Gaussian upper and lower bound for the fundamental solutions of (jl.lj) . The 
third theorem is an application of the second one in the more traditional setting of L p,q 
conditions on the potential. It will show that our conditions are genuinely much broader 
than the traditional ones. 

It should be made clear that there is no claim on uniqueness in any of the theorems. In 
the absence of uniqueness how does one define the fundamental solution? This is possible 
due to the uniqueness of problem (jl.lj) when the potential V is truncated from above. This 
fact is proved in Proposition 2.1 below. Consequently we can state 

Definition 2.2. The fundamental solution Gy is defined as the pointwise limit of the 
(increasing) sequence of the fundamental solution Gy i where Vi = min{V, i} withi = 1, 2, .... 

We remark that Gy thus defined may be infinity somewhere or everywhere. However we 
will show that they have better behavior or even Gaussian bounds under further conditions. 

Theorem 2.1. (i). Suppose M-l)) with some uq>Q has a positive solution. Then 

V = Af-\Vf\ 2 -d t f, 

with e~ f E L 2 (R n x R+). Moreover f E L} oc (R n x R+) if\nu E L} oc (R n ). 
(ii). Suppose 

V = Af-\Wf\ 2 -d t f 

for some f such that e~l E L 2 (R n x R + ). Then the equation in M.l\) has a positive L 2 
solution for some Uq E L 2 (R ri ). 

Theorem 2.2. (i). Suppose 

V = Af-a\Vf\ 2 -d t f 

for one given a > 1 and f E L°°(R" x R + ). Then Gy has Gaussian upper bound in all 
space time. 

(ii). Under the same assumption as in (i), if gi/i* |V/| 2 E L°°(R n x R + ), then Gy has 
Gaussian lower bound in all space time. 

(Hi). Under the same assumption as in (i), suppose Gy has Gaussian lower bound in all 
space time. Then there exists h > such that 

2 b *|V/| 2 eL 00 (R"xR+). 

Remark 2.1. At the first glance, Theorem 2.1 may seem like a restatement of existence 
of positive solutions without much work. However Theorem 2.2 shows that if one just puts 
a little more restriction on the potential V, then the fundamental solution actually has a 
Gaussian upper bound. Under an additional but necessary assumption, a Gaussian lower 
bound also holds. Even the widely studied potential a/|x| 2 in R™ can be recast in the form 
of Theorem 1.1, as indicated in the following 
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Example 2.1. For a real number b, we write / = felnr with r = \x\. Then direct 
calculation shows, for r ^ 0, 

*^ = 4/-|v/|> 

Let a = &(n — 2 — 6). Then it is clear that the range of a is (—oo, (n — 2) 2 /4]. In this 
interval (1.1) with V = a/\x\ 2 permits positive solutions. This recovers the existence part 
in the classical result [BG]. Highly singular, time dependent examples can be constructed 
by taking / = sin(^^) e.g. 

Moreover the corollaries below relate our class of potentials with the traditional "form 
bounded" or domination class ()2.1|) . (see also jHIj below). In the difficult time dependent 
case, Corollary 1 shows that potentials permitting positive solutions, can be written as the 
sum of one form bounded potentials and the time derivative of a function almost bounded 
from above by a constant. 

Remark 2.2. From the proof, it will be clear that under the assumption of part (i) of 
Theorem 2.2, one has 



/ 



G a y(x, t] y, s)dy < C < oo. 



This is one of the main assumptions used by Sturm [St] in the time independent case 
(Theorem 4.12). If a = 1, then the conclusion of Theorem 2.2 may not hold even for time 
independent potentials. See [St]. Also note that this theorem provides a nearly necessary 
and sufficient condition on certain sign changing potential such that the corresponding 
heat kernel has Gaussian upper and lower bound. The only "gap" in the condition is the 
difference in the parameters of the kernels #1/4 and g^. It is well know and easy to check if 
|V/| G L p ' q with - + ~ < 1 and / = outside a compact set, then g>&* |V/| 2 is a bounded 
function for all b > 0. 

Corollary 1. Let V G L 2 (R n x (0,oo)). 

(a) . Suppose 

(2.1) f [ V<?< f [ |V0| 2 + 6 f 

Jo J Jo J Jo 

for all smooth, compactly supported function <fi in R n x (0, T) and some b > and T > 0. 
Then M.l)) has a positive solution when u > and moreover 

V = Af-\Vf\ 2 -dJ. 

(b) . Suppose V = Af — |V/| 2 then V is form bounded, i.e. it satisfies K2.1)) . 

(c) . Let V be a L 2 potential permitting positive L 2 solutions for (1.1), then V can be 
written as the sum of one form bounded potentials and the time derivative of a function 
almost bounded from above by a constant. 



In the next corollary, we consider only time- independent, nonnegative potentials. Here 
the definition of V = Af — |V/| 2 is slightly different from that of Definition (1.1) since 
we do not have to worry about time derivatives. One interesting consequence is that these 
class of potentials is exactly the usual form boundedness potentials. 
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Corollary 2. Suppose <V £ L 1 (R n ). Then the following statements are equivalent. 

(1) . For some f G L} oc (H n ) and a constant b > 0, 

V = Af-\Vf\ 2 + b. 

This mean there exist Vj G L°° such that Vj — > V in L 2 (R 1 ) and Vj = Afj — |V/j| 2 + b for 
some fj G W 2 ' 2 (R n ), j = 1,2,.... 

(2) . 

J V<p 2 dx < J \V(p\ 2 dx + bj <p 2 dx. 
for all smooth, compactly supported function (f> in R n and some b > 0. 

Remark 2.3. Condition (2) in the above corollary just means that the bottom of the 
spectrum for the operator —A — V is finite. This condition is the same as those given in 
[CM] and [GZ]. 

It is a fact that most people feel more familiar with the case when the potential V is 
written as L p,q functions. Also there may be some inconvenience about the presence of the 
nonlinear term |V/| 2 in the potential in Theorem 2.2. Therefore in our next theorem, we 
will use only L p,q conditions on / without nonlinear terms. 

Theorem 2.3. Suppose 

(a) f G L°°(R n x R+); 

(b) f = outside a cylinder B(0, R ) x [0,T ], Rq,T > 0; 

(c) IV/ 1 G L p ' q (K n x R+) with 2 + 2 < i; 

(d) V = Af -d t f G L 2 (R n x R+). q 

Then there exists a constant Aq depending only on n, p, q such that the following state- 
ments hold, provided that 

I |V/| || LP'i (R»xR+) < A . 

(i) The kernel Gy has Gaussian upper and lower bound in all space time. 

(ii) The kernel Gg t f has Gaussian upper and lower bound in all space time. 

Remark 2.4- If V is independent of time, then Theorem 2.3. reduces to the well known 
classical fact: 

if a potential V is the derivative a of a small L n+e function, then Gy has Gaussian upper 
and lower bound. (See [LS]) e.g. 

In the time dependent case our result is genuinely new due to the presence of the term d t f. 
Let us mention that some smallness condition on the potential is needed for the existence 
of Gaussian bounds for Gy. This is the case even for time independent, smooth potentials 
due to the possible presence of ground state. 



2.2. Preliminaries. In order to prove the theorems we need to prove a proposition con- 
cerning the existence, uniqueness and maximum principle for solutions of under the 
assumptions that V is bounded from above by a constant. The result is standard if one 
assumes that the gradient of solutions are L 2 . However we only assume that solutions are 
L 2 . Therefore a little extra work is needed. 
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Proposition 2.1. Suppose that V G L 2 (R™ x R + ) and that V < b for a positive constant 
b. Then the following conclusions hold. 

(i) . The only L 2 solution to the problem 

Au + Vu-d t u = 0, R n x (0, T), T > 0, 
u(x, 0) = 0, 

is zero. 

(ii) . Let u be a solution to the problem in (i) such that u(-,t) G L 1 (R n ) ; u G L 1 (R" x 
(0,T)) and Vu G L^R" x (0,T)). T/jen u is identically zero. 

(Hi). Under the same assumptions as in (i), the following problem has a unique L 2 
nonnegative solution. 

' Au + Vu-d t u = 0, R n x (0, T), T > 0, 
u(.,0) = uo(0 >0, u eL 2 (R n 

(iv). Under the same assumptions as in (i), letu be a L 2 solution to the following problem 

Au + Vu- d t u = f, R n x (0, T), T > 0, 
u(-,0) = 0. 

Here f < and f e L 1 (R n x (0,T)). T/ien u > m R n x (0,T). 

Proof of (i). Let u be a L 2 solution to the problem in (i). Choose a standard mollifier p 
and define, for j = 1, 2, 

Uj(x,t)=j n J p(j(x-y))u(y,t)dy = J pj(x - y)u(y,t)dy. 

Then Vuj and A-u^ exist in the classical sense. From the equation on u, it holds 

Auj + J pj(x - y)V(y,t)u(y,t)dy - d t Uj = 0, 

where d t Uj is understood in the weak sense. 
Given e > 0, we define 

hj = \Ju 2 + e. 

Let < <fi G C^°(R ri ). Then direct calculation shows 
J (fihj(x, s)\ t dx = 



a 



u -Aw • /"* /" u (x s) f 
3 3 (x, t)(p(x)dxds + / / ^ = — / pj(x — y)V(y,t)u(y,t)dydxdt 

/u 2 + e Jo J Jtf + t J 



= n+T 2 . 
Using integration by parts, we deduce 



Ti = — I I ] ^ Uj][ s)dxds + [ [ y==4>{ x , s)dxds 



Uj + e Jo J j ' J«j + e 



2 _|_ e ./n ./ "7 + r /,,2 



u 2 + e 
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Since the sum of the first two terms on the righthand side of the above inequality is non- 
positive, we have 

\Vuj\ |V(/>| 



Ti < 



u. 



Taking e to zero, we obtain 



-dxds. 



< 



Vuj\ | V(j)\dxds + / / 4>(x)\ / pjV + u(y, s)dy\dxds 



\Uj 



4>(x)— —Ax, s) I pj(x — y)V u(y, s)dydxds. 



Here and later we set s) = if Uj(x, s) = 0. 

Next, since u,V e L 2 (R" x (0,T)), one has uV G L^R" x (0,T)). From the equation 

Au + Vu- d t u = 

one deduces 

u(x,t) = / / G (x,t;y,s)(Vu)(y,s)dyds. 



Here, as always, Go is the fundamental solution of the free heat equation. Hence u(-,t) G 
L 1 (R") and u G L 1 (R™ x (0,T)). Therefore, for any fixed j, there holds 

|Vuj| G L x (R n x (0,T)). 

Now, for each R > 0, we choose so that <fi = 1 in B(0,R), = in B(0,R + l) c and 
|V0| < 2. Observing 

J j |V0|cfeds -> 0, R->oo, 

we deduce, by letting i? — > oo, 
(2-2) 

y < y J \J PjV +u (Ui s )dy\dxds— J J j-^(x,s) J pj(x—y)V~u(y, s)dydxds. 

By the fact that V~u, V + u G L 1 (R n x (0,T)), we know that 

- y)V~u(y, -)dy -> V~u(-, •), 



in L 1 (R n x (0,T)). Since y^yO^, s) is bounded and converges to ^(a;, s) a.e. in the support 



of u, we have 



M 



< 



u- f V~u 2 

[y^-Axis) j pj(x - y)V~u(y,s)dy — — ]dxds 

o J \ u j\ J \ u \ 

1 [ (1L-U, v _ udxds 
o J \ u j\ \ u 







U j 



U : 



(pj *V u — V ujdxds 



0. 
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Substituting this to (|2.2|) we deduce, by taking j — ► oo, 

/"* f f* f V~u 2 

\u(x,t)\dx < J J V + \u(x, s)\dxds — J J — j — — (x,s)dxds. 

Therefore 

\u(x, t)\dx < J J \u(x, s)\dxds \\V + \\ 00 . 

By Grownwall's inequality u(x,t) = a.e. This proves part (i). 
Proof of (ii). 

Notice that the only place we have used the L? boundedness of u is to ensure that 
Vu G L 1 (R n x (0, T)). But this a part of the assumptions in (ii). Therefore (ii) is also 
proven. 

Proof of (in). The uniqueness is an immediate consequence of part (i). So we only need 
to prove existence. This follows from a standard limiting process. For completeness we 
sketch the proof. 

Given k — 1, 2, .. let be the truncated potential 

Vk = sup{V(x,t), -k}. 

Since V%. is a bounded function there exists a unique, nonnegative solution u k to the following 
problem. 

Au k + V k u k - d t u k = 0, in R n x (0, oo) 
u(x, 0) = u (x) > 0, x G R n , u G L 2 (R n ). 

By the standard maximum principle, {u k } is a nonincreasing sequence and 

- J u 2 k \^dx+ J J \Vu k \ 2 dxdt = J J v k u\dxdt < J J V k + u\dxdt < \\V + \\ 00 J J u 2 k dxdt. 
By Grownwall's lemma, we have 

»T r r r r 



J J | Vu k \ 2 dxdt + J ul(x,T)dx < J ul(x)dx + ||t' + ||oo J Uq\ 



x)dx e 



It follows that u k converges pointwise to a function u which also satisfies the above inequal- 
ity. Let be a test function with compact support. Then 



"T r rT r rT 



(u k (J))\ dx — / / u k (ptdxdt — / / u k A(f)dxdt — / / V k u k (j)dxdt = 
Jo J Jo J Jo J 

Since {u k } is a monotone sequence and also since |VfcMfc| < G L^R" x (0,T)), the 

dominated convergence theorem implies that 

/{ucfy^dx — I I u(f>tdxdt — / / uA(f)dxdt — / / Vucf)dxdt = 0. 
Jo J Jo J Jo J 

This shows that u is a nonnegative solution. It is clear that u is not identically zero since 
uq is not. This proves part (iii) of the proposition. 

Proof of Part (iv). 
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Let V k be a truncated potential as in part (iii). Since V k is bounded, the standard 
maximum principle shows that there exists a unique, nonnegative solution to the following 
problem. 

' Au k + V k u k - d t u k = f < 0, R" x (0, T), T > 0, 
« fc (-,0) = 0. 

Moreover {u k } forms a decreasing sequence. Since V k is a bounded function, the standard 
parabolic theory shows that 



Therefore 



/u k (x,t)dx = / / V k + u k dxds — / V k u k dxds + / / fdxds. 
Jo J Jo J Jo J 

J u k (x,t)dx < \\V + \\oo J J Ukdxds + J J fdxds. 

This implies 

J u k (x,t)dx + J J u k dxds KC^WV+W^WflU). 

It follows that 

J u k (x,t)dx + J J u k dxds + J j V k -u k dxds < C(t, ||^ + ||oo, H/lli). 
Let iu be the pointwise limit of the decreasing sequence {u k }. Then we have 
J w(x,t)dx + J J wdxdsKC^WV+lU ll/HO. 

It is straight forward to check that w is a nonnegative solution to the problem 

Aw + Vw-d t w = f<0, R n x (0, T), T > 0, 
w(-,0) = 0. 

Hence 

|A(w-H) + y(w-M)-9 t (w-tt)=0, R n x (0,T), T>0, 
\ («;-«)(., 0) = 0. 

Recall that -u is assumed to be a L 2 solution and that V £ L 2 . We have that Vu G L 1 
and consequently G L 1 (R n ) and u G L 1 (R n x (0,T)). Now by Part (ii) of the 

proposition, we deduce w = u since w is also L 1 . Hence u > 0. This finishes the proof of 
the proposition. □ 



2.3. Proofs of Theorems. 



Proof of Theorem 2.1 (i). 
For j = 1,2..., let 

Vj = min{V(x,t),j}. 

Since Vj is bounded from above, Proposition 2.1 shows that there exists a unique solution 
Uj to the following problem. 

Auj + VjUj - d t Uj = 0, (x, t) G R™ x (0, oo) 
Uj(x, 0) = u (x), x G R ra . 
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Notice that Uj — Uj_i is a solution to the problem 

A{u,j — Uj-i) + Vj(uj — Uj-i) — d t (uj — Uj-i) = (Vj-i — Vj)uj-i, (x, t) G R n x (0, oo) 
(uj -uj_i)(a:,0) = 0, x G R™. 

Notice that 

(tj-i - ^K-i < 0, (^-i - ^K-i e L^R" x (0, T)), T > 0. 
We can then apply Proposition 2.1 (iv) to conclude that 

Uj > Uj-i. 

Moreover 

A(u - uj) + Vj(u - Uj) -d t (u- Uj) = (Vj - V)u, (x, t) G R™ x (0, oo) 
{u-Uj)(x,0) = 0, i6R B , 

with 

(Vj-V)u<0, (Vj — V)u G L 1 (R n x (0,T)), T > 0. 
By Proposition 2.1 (iv) again we know that 

U > Uj. 

Therefore {uj} is a non- decreasing sequence of nonnegative functions bounded from above 
by a L 2 function. Let w be the pointwise limit of Uj. The w is L 2 and \VjUj\ < \Vu\ G 
L 1 (R n x (0,T)), T > 0. By the dominated convergence theorem, it is straight forward to 
check that w is a nonnegative L 2 solution to the equation 

Aw + Vw- d t w = 0, (x, t) G R n x (0, oo) 
w(x,0) = uq(x), x G R n . 

Fixing j, for any fc = 1, 2, Let 

Vjk = max{Vj(x, t), —k}. 

Since Vjk is bounded, the following problem has a unique L 2 solution. 

Aujk + VjkUj k - d t Ujk = 0, (x, t) G R n x (0, oo) 
Uj(x, 0) = m (x), x G R n . 

Due to the fact that {Vjk} is a decreasing sequence of k, the maximum principle shows 
that {uj k } is also a decreasing sequence of k. Since Vj k u jk G L 2 (R n x (0,T)), T > 0, the 
parabolic version of the Calderon-Zygmond theory shows 

Au jk , d t u jk G L 2 (R n x (0, T)),T> 0. 

Since 

< Ujk — Uj < Uji — Uj, k — 1, 2, 3, 
< ty — < w — u±, 
we can apply the dominated convergence theorem to conclude that 

lim / / (v,jk — Ujfdxdt = 0, lim / / (w — Uj) 2 dxdt = 0. 
k ^°°Jo J Jo J 
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Therefore we can extract a subsequence {uj^} such that 

lim / (w — Ujk) 2 dxdt = 0. 
Jo J 3 

Hence there exists a subsequence, still called {u/j^.} such that 

- > w a.e. 

Recall that Uq > 0, u ^ and is bounded. It is clear that Ujk j > when t > 0. 
Now we define 

fj = -\nu jkj , f = -\ogw. 

Then 

\},. a./; ; - Yfr-0,1). 

Clearly fj — > f a.e. and Vjk- — > V in L 2 as j — > oo. By Definition 2.1, this means 

V = Af-\Vf\ 2 -d t f. 
It is clear that e~^b = w is L 2 by construction. 
Proof of Theorem 2.1 (ii). 

By assumption, there exist sequences of functions {Vj} and {fj} such that 
\\Vj\\ L 2 <C, |^-|<l^+i|, WVj-VWv^oa, 
Vj = A/ 3 - - |V/j| 2 - d t fj. 

Then for w,,- = e~h , we have 

A«j + VjUj — dtUj = 0. 
We will show that ||mj||l2 is uniformly bounded. To this end, we observe that 

A(uj - u j+ i) + Vj(uj - Uj+i) - d t {uj - u j+ i) = -(Vj - V j+1 )u j+1 > 0. 

Recall from Definition 2.1 that fj(x,0) is independent of j. Hence Uj(x, 0) = Uj + %(x, 0). 
Therefore < Uj < Uj+i. By the assumption that fj — *> f a.e., we know that Uj = e~h — > 
e~l a.e. Note that e~? G L 2 (R™ x (0, oo)). Hence ||Wj||L a is uniformly bounded. 

By weak compactness in L 2 , there exists a subsequence, still called {uj} such that Uj 
converges weakly to a L 2 function which we will call u. Observe that, for any compactly 
supported test function 0, there holds 

< ||^(^--V)0|| L i + ||K-- M )\/0|| L i 

< ll%IU 2 ll^--^IU 2 1101k- + \\{uj-u)V(j)\\ L i. 

Hence 

\\ujVj(p-uV(p\\ L i -> 

when j — > oo. From here it is easy to check that u is a nonnegative solution to (jl.lj) with 
Uq = e~^^ x ' ' as the initial value. Note the u$ is independent of j. If u G L 2 (R n ) then we 
are done. Otherwise, we can selection a L 2 function dominated by uo to serve as the initial 
value. □ 

Next we will provide a 
Proof of Theorem 2.2 (i). 

We will use an idea based on an argument in [St] where the heat equation with some 
singular, time independent potentials are studied. 



LOCAL REGULARITY 13 

By virtue of Proposition 2.1, the fundamental solution Gy is defined as the limit of 
fundamental solutions of the equation in (jl.lj) where V is replaced by nonsingular potentials. 
Therefore we can and will assume that V is smooth in this subsection. The constants 
involved will be independent of the smoothness. 

Since, by assumption 

(2.3) V = Af-a\Vf\ 2 -d t f, 
one has 

aV = A(af)-\V(af)\ 2 -d t (af). 
Writing F = e~ a f , it is easy to show that 

(2.4) AF + aVF - d t F = 0. 

Let us denote the fundamental solution of the equation in ()2.4j) by G a y. Since / is bounded, 
we know that F is bounded between two positive constants. Therefore it is clear that 

inf F f ^ . , sup F 

2.5 < < / G aV {x, t; y, s)dy < -f--, 

sup r J mi t 

for all x G R" and t > s. Here inf F and sup F are taken over the whole domain of F. 
By Feynman-Kac formula and Holder's inequality, for a given G C^°(R n ), there holds 

G v {x,t V y,s)(j ) { y )d y \ < [ J G aV {x,t- y ,s)d y ] 1/a [ j G (x, t; y, s)\<f>{y)\ a l^dy\ (Q ~ 1)/a 

By (12.5)) . we deduce 

l lot 

cs 

G v {x,t;y,s)<f>{ y )dy\ < - ^ ( °_ 1)w/(2a) H\\a/{a-i), 
where So = fr^f?- The norm on means the L a ^ a-1 ^(R n ) norm. Hence 

l/a 

(2-6) ||G ! y(-,t;-,s)|| a/(a _ 1))00 < ^_^ ( °_i )n/(2a) - 

Here and later the norm || ■ || Pj g stands for the operator norm from L p (R n ) to L 9 (R n ) for 
p, q between 1 and oo. 

Without loss of generality we assume that a /(a — 1) is an integer. This is so because oth- 
erwise we can choose one a\ G (1, a) such that a.i/{a\ — 1) is an integer. Then interpolating 
between G a y and Go by Feynman-Kac formula again, we know that 

G aiV (x, t; y, s)dy < C(F, a, a x ). 

Then we can just work with G ai y instead of G a y in the above. 
Using the reproducing property of Gy we deduce 

(2.7) ||G y (-,t;-, S )|| 1 , 00 <n j m =1 ||G' y (-, S + ^; •, Vi)|| Pj , gj , 

where 

m = a/(a — l), Pj — m/(m — I + 1), qj—m/(m — l), tj — s + ((i — s)j/m). 
For each j between 1 and m, we apply the Riesz-Thorin interpolation theorem to deduce 



!- a j \\n I 4- . 4- Ml A 



\Gy(-,tj; -j^-ijllpj,^ < \\Gy(-,s + tj; •,^-i)|li m /( m _i) \\Gv{-,tj] -,£/-i 



3 

I m,oo ' 
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Here the parameters are determined by the following relations 

1 1 — A,- Xj 1 1 — Xj Xj 
Pj 1 m ' qj m/(m — 1) oo ' 

A, = 1^1. 
m — 1 

It follows that 

(2-8) \\G v (-,tj] ■,tj-i)\\ Pjtqj < \\Gv(-,tj] -, tj-i)|| m ,cx3- 

Substituting (|2.6|) to (|2.8|) . we deduce, after noticing that t,- — tj-i = (t — s)/m, 

cs 1 /« m (a-l)n/(2a) 



|Gy(-,t i ; ^tj^Wp^q. < _ s \( a _i) n /(2a) " 



This and ()2.7|) imply that 



|GVM;^)lkoo<^ 



/ a r? /2 

cs n m ' 



(t - s) n l 2 ' 

Here we just used the relation m = a/ (a — 1). This yields the on-diagonal upper bound 

csn /[a ~ 1] (a/(a - l)) n/2 

(2.9) G v (x, t; y, s) < ^ ( /^ /2 ^ ■ 

In order to obtain the full Gaussian bound, we observe that, for any p > 1, the Feynman-Kac 
formula implies 

(2.10) G v (x, t; y, s) < [G pV (x, t; y, s)] 1/p [G (x, t; y, s)] iP ~ 1)/p . 
Notice also 

pV = p(Af - a\Vff - d t f) = A(pf) - -\V(pf)\ 2 - d t (pf). 

p 

Taking p = (1 + a)/2, then - > 1. Therefore pV also satisfies the condition of Theorem 
2.2 (i). Hence, the on-diagonal bound ()2.9|) holds for G v y. i.e., there exists a constant 
C(a,e sup/ - inf/ ) such that 

„ , . C(a,e su P^ inf/ ) 

G pV {x,t;y,s) < ^ _ . 

Substituting this to the inequality ()2.10|) . we obtain the desired Gaussian upper bound for 

Gy. 

Proof of (ii). 

In this part we prove the Gaussian lower bound. We will follow Nash's original idea. The 
novelty is a way of handling the potential term even if it is very singular. The main idea is 
to exploit the structure of the potential when it is written as a combination of derivatives. 

Since the setting of our problem is invariant under the scaling, for r > 0, 

V r (x,t) = r 2 V(rx,r 2 t), f r (x,t) = f(rx,r 2 t), u r (x,t) = r 2 u(rx,r 2 t), 

we can just prove the lower bound for t = 1 and s = 0. We divide the proof into three 
steps. 

Step 1. Fixing x G R n , let us set 

u(y,s) = G v (y,s;x } 0), 
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H(s) = J e- nlyl2 \nu(y,s)dy. 
Differentiating H(s), one obtains 



H'(s) = I e-^ 9 -4^4 dy 
u(y,s) 



- J V(^—^-)Vudy + J e'^ y \ 2 V(y,s)dy. 



Estimating the first term on the righthand side of the above inequality as in [FS], section 
2, one arrives at 

(2.11) H\s) > -C + \ J e~ 7rlyl2 \V\nu(y,s)\ 2 dy + J e^Vfe/, s)dy. 

Here C is a positive constant. Since, 

V = Af-a\Vf\ 2 -d t f 1 

we know that 

f(x,t) = J G (x,t;y,0)f(y,0) 

-J J G (x,t;y,s)V(y,s)dyds- a j J G {x,t;y, s)|V f\ 2 dyds. 
By our assumption 

J J G (x,t; y, s)|V f\ 2 dyds < oo. 
Hence the boundedness of / implies that 

m(x,t) = — / / G (x, t; y, s)V(y, s)dyds e L°°. 



o 



Moreover 

(2.12) Am-d t m = V. 
Therefore 

J e-^ y? V{y,s)dy = J e^ |y|2 [Am - d s m)(y, s)dy 
= J [Ae-^mdy - d s J e-^ 2 m{y, s)dy 

>-C-d s J e~ 7rlyl2 m(y,s)dy. 

Here we have used the boundedness of m. Substituting the above to the righthand side of 
(|2.1ip . we obtain 

(2.13) H'(s) >-C + ^ J e- wlyl2 \V\nu(y,s)\ 2 dy- M'(s), 
where 

(2.14) M(s) = j e- nlyl2 m(y,s)dy. 
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Step 2. By Poincare's inequality with e _7r ' y ' 2 as weight, we deduce, for some B > 0, 

H'{s)>-C + B J e-^ y? [\nu{y,s)-H{s)] 2 dy-M'{s). 

Next, observe that (lnw — H(s)) 2 /u is non-increasing as a function of u when u is between 
e 2+H( s ) an( j ^ ^j go f rom fae Gaussian upper bound, 



sup u(y, s) < K < oo. 

1/2<s<1 



Therefore 



(2.15) H'(s)>-C + CBK-\\nK - H{s)f [ e~ n ^ 2 u(y, s)dy — M'(s). 

J u(y,s)>exp(2+H(s)) 

Using the Gaussian upper bound again, we know that H(s) < C for some C > and 
that 

f e- nlyl2 u(y,s)dy> [ e^ 2 u(y, s)dy - ce 2+H{s) 

J u(y,s)>exp(2+H(s)) J 

(2.16) >e~ m 2 [ u(y,s)dy-ce 2+H ^ 

J\v\<r 



-irr 



\y\<r 

y, s)dy - I u(y, s)dy] - ce 2+H{s) . 

J\y\<r 

We aim to find a lower bound for the righthand side of ()2.16|) . By ()2.12j) . 

V = Am — dtm > Am — |Vm| 2 — dtvn = V\. 

Write h = e~ m . Then 

Ah + Vih - d t u = 0. 

Since m e we know that h is bounded between two positive constants. Observe that 



h(x,t) = J G Vl (x,t;y,s)h(y,s)dy. 

Hence 

< ci < J G Vl (x, t; y, s)dy < c 2 . 
By the maximum principle, we have 

(2.17) J G v (x,t;y,s)dy > J G Vl (x, t; y, s)dy > c x > 0. 

Recall that u(y, s) = Gy(y, s; x, 0). Substituting (|2.17j) to (|2.16|) and applying the Gauss- 
ian upper bound on u(y, s), we deduce 

(2.18) f e-^ 2 u{y, s)dy > ce^d - ce 2+H{s \ 

J u{y,s)>exp{2+H{s)) 

when r is sufficiently large. Substituting ()2.18|) to (|2.15)L we arrive at 

(2.19) H'(s) >-C + CBK-\\nK- #(s))V" r2 c x - ce 2+H{s) - M'(s). 

We claim that H(l) > —cq for some sufficiently large Co > 0. Suppose otherwise, i.e. 
H(l) < -c . From (l2~THj) . for some C > 0, 

H'{s) > -C-M\s). 
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Hence 

HQ) - H(s) > -C(l - s) - (M(l) - M(s)). 

Therefore 

H(s) < H(l) + C(l - s) + (M(l) - M(s)) < -c /2, 
when Co is chosen sufficiently large. It follows from ()2.19jl that 

H'(s) > ~c 1 + c 2 H 2 (s)-M'(s). 

This shows 

(H(s) + M{s))' > -c 3 + c A {H{s) + M{s)f. 
From here, one immediately deduces 

H(l) >-A, A > 0. 

The claim is proven. Thus 



where \x — y\ < 1. 

Using the reproducing property of Gy and Jensen's inequality, we have, when \x — y\ < 1, 
In Gy(x, 2; y, 0) = In J Gy{x, 2; z, l)Gy(z, 1; y, 0)dz 

> J e- 7Tlvl2 \nG v (x,2;z,l)dz + J lnG v (2, 1; y, 0)dz > -C. 

This proves the on-diagonal lower bound. The full Gaussian lower bound now follows from 
the standard argument in [FS]. □ 

Proof of (in). 
Since, 

V = Af-a\Vff-d t f, 

we have 

V+(a-l)|V/| 2 = A/-|V/| 2 -cV. 
Let u = e~f , by direct calculation, 

Au + Vu- d t u + (a - 1)|V/| 2 m = 0. 

Hence 

u{x,t) = / G v (x,t;y,0)u(x,0)dy+ {a - 1) / / G v (x, t; y, s)\Vf\ 2 {y, s)u{y, s)dyds. 



Since / G L°°, we know that u is bounded between two positive constants. If, by assump 
tion, Gy has a Gaussian lower bound, then, for some b > 0, we have 

sup u 

g b (x,t;y, s)\V f\ 2 (y, s)dyds < C^— 

mi u 

This completes the proof of part (iii) of Theorem 2.2. □ 

Proof of Theorem 2. 3. 
(i). We write 

y 1 = 2(A/-^/-3|V/| 2 ), ^2 = 6|V/| 2 . 
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0. Since 



1 /2 

G v (x,t;y,s) < [G Vl {x,t;y,s)\ [Gy 2 {x,t;y,s 



>] 



1/2 



Observe that 



\/ 1 = A(2/)-a t (2/)-^|V(2/)| 2 . 



Hence, by Theorem 2.2, we know that Gy 1 has Gaussian upper bound. Under the smallness 
assumption on the L p,q norm of |V/| 2 in the theorem, it is well known that Gy 2 also has a 
Gaussian upper bound. Therefore Gy has Gaussian upper bound. 
In order to prove the Gaussian lower bound, we observe that 



Under our assumption on the L p ' q norm of |V/|, it is straight forward to check that 



Hence Theorem 2.2 (ii) shows Gy 3 has Gaussian lower bound. Clearly this Gaussian lower 
bound of Gy 3 is also a Gaussian lower bound of Gy by the maximum principle. This proves 
part (a). 

(ii). Clearly we can choose Aq sufficiently small so that all the following kernels have 
global Gaussian upper and lower bound: 

(2.20) G2V1 Gy/2, G2A/, Ga//2- 

The bounds on the first two kernels follow from part (i). The bounds on the last two kernels 
follow from standard theory since Af = div(Vf) with V/ has a small norm in the suitable 
L p,q class, (see |LSj e.g.) 
Now observe that 



We close this section by giving proofs of the corollaries. 
Proof of Corollary 1. 

(a). Let Vk = min{V(x, t), k}, k = 1,2,.... Then (II. lj) with V replaced by V& has a 
unique solution. 

Let J(t) = J nn ul(x,t)dx. Then 



V = Af - d t f >V 3 = Af - d t f - 2|V/| 2 . 



gi/**\Vf\ 2 eL 



-dtf — Af - d t f — Af — V - Af, 
V = Al_dJ 
2 2 2 



By Feynman-Kac formula 
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By our assumption on V, 

J'(t) < 2bJ(t) 

which implies 

/ u 2 k (x,t)dx < / ul(x)dxe 2bt . 
Jr.™ Jn n 

Therefore if uq G L 2 (D), we conclude that Uk(x,t) increases to a finite positive limit 
u(x,t) as k — > oo, for all t and for a.e. x. Moreover u(-,t) G L 2 (R n ). We show that the 
above u is a positive L 2 solution to 

Since Uk is a solution to (jl.lj) with V replaced by Vk, for any ip G C^°(R n x (0,T)), we 
have 

rt2 r rti r rti r 

(ukip)\lldx — / / Ukiptdxdt — / / UkAipdxdt — / / VkUk^dxdt = 

y j<i y j<i y 

for all ti,t 2 G (5, to)- 

By our assumption | V^Wfe | < \Vu\ G L 1 (R n x (0,T)). Taking A; — > oo and using the 
dominated convergence theorem, we obtain 

[uip)\lldx — / / uvp t dxdt — / / uAipdxdt — / / Vuipdxdt = 0. 

Ai y y<i y y*i y 

This shows that u is a positive solution to (jl.lj) . By Theorem 2.1 

V = Af-\Wf\ 2 -d t f. 

(b) . Suppose 

^ = A/-|V/| 2 . 

Due to the L 2 convergence, it suffices to prove that Vj in Definition 2.1 satisfies (j2.1j) . Let 
be a test function, then 

J V j( f) 2 dxdt = J J [Afj - \Vfj\ 2 ](j) 2 dxdt 

/oo r poo r 

J V/jV0 0dxdt - J J \Vfj\ 2 (j) 2 dxdt 

/OO /* 
J \V<p\ 2 dxdt<p. 

(c) . The statement is self-evident by part (b) and Theorem 2.1. □ 
Proof of Corollary 2. 

Suppose V = Af — | V/| 2 + b. Then, by the same limiting argument as above, we have 
Vtfdx = lim / [Afj - \VfjWdx + b [ <p 2 dx 

j^OD J J 



lim ( - 2 J VfjVcj) <pdx - j |V f '<p 2 'dx \ +b j <p 2 dx. 



Therefore 

Also by part (a) of Corollary 1, (jl.lj) has a positive solution when uq > 0. 



20 



QI S. ZHANG 



On the other hand, suppose V satisfies 

j V4> 2 dx< J \V<p\ 2 dx + b J <fdx. 
Write Vj = min{V, j} with j = 1, 2, .... Then 

{Vj-b)4) 2 dx < [ \V(j)\ 2 dx. 



Notice that Vj — b is a bounded function. Hence we can apply Theorem C.8.1 in [51] to 

;uch that 

A Uj + {V, - b) Uj = 0. 

V j = Af j -\Vf j \ 2 + b. 
V = Af-\Vf\ 2 + b. 



conclude that there exists Uj > such that 



Writing fj = — InUj, we have 
By definition, this means 



□ 



3. HEAT BOUNDED FUNCTIONS AND THE HEAT EQUATION 



Here we introduce another class of singular functions that has its origin in the Kato type 
class. As mentioned in the introduction, a function is in a Kato type class if the convolution 
of the absolute value of the function and the fundamental solution of Laplace or the heat 
equation is bounded. Here we generalize this notion by a simple but key stroke, i.e., we 
delete the absolute value sign on the function in the definition of the Kato class. More 
precisely, we have 

Definition 3.1. Let / = f(x, t) be a local L l function in space time and Go be the standard 
Gaussian in R n . We say that / is heat bounded in a domain Q C R" x R 1 if 

G *f(x,t)= / G (x,t;y,s)f(y,s)dyds 

JO JR, n 

is a bounded function in Q. 

We say that / is almost heat bounded in a domain Q C R n x R 1 if 

G *f(x,t) eL"(n) 

for all p > 1. 

Example. The function V(x) = a x ^ 2 1) is not heat bounded but is almost heat bounded 
in R". Here a is a nonzero constant. 

In the next two propositions, we provide a comparison between the heat bounded class 
and more familiar classes of functions. 

Proposition 3.1. Suppose, in the distribution sense, V = dfjf with f 6 n p >iL p (R ?1 x 
(0,T)). Then V is almost heat bounded in R n x (0, T). 

Proof. 
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Let G be the free heat kernel in R™ x (0, oo). By the assumption on /, the function 
u = u(x,t), defined by 

u(x,t)= / G (x,t;y)f(y,s)dyds 
Jo Jr" 

is a solution to the equation 

Au(x, t) - u t (x, t) = -f(x, t), x G R n , t > 0; 
u(x,0) = 0. 

By the parabolic version of the Calderon-Zygmond inequality (see [Lieb] e.g., we know that 

u e W 2 ' p {R n ), Yp > l. 

Hence 

d%\ [ G (x,t;y)f(y,s)dydseL^R n ), Vp > 1. □ 
Jo Jr™ 

Proposition 3.2. Suppose, < V G L] oc is form bounded in D x [0, T\. i.e. 

f [ V^Kb, f [ |V0| 2 + 6 2 f [ 2 . 
Jo Jd Jo Jd Jo Jd 

for all smooth, compactly supported function <fi G D x [0, T] C R n x [0, T] . Then V is almost 
heat bounded in D x [0, T]. 

Proof 

We will only consider the case when D = R n . The other cases follow from the full space 
case by a standard comparison method. 

Since one can consider cV with c sufficiently small otherwise, we can choose the constant 
bi in the definition of form boundedness to be 1/2, i.e. we assume that 

rT r -| pT p pT 



Jo Jd z Jo Jd Jo Jd 



for all smooth, compactly supported function 0. 
Let Uk be the solution of 



(3.1) 



Au k + V k u k - (u fc )t = 0, in R" x (0, oo), V G Lf oc (R n x (0, oo)) 
u k (x, 0) = u (x) > 0, x G R n , M G L 2 (R n ). 



Here V k is the truncated potential V k = min{V, k} with k being positive integers. Clearly 
V k < Vfe+i. 

We show that u k converge pointwise to a locally integrable function. 
Let J(t) = f D ul(x,t)dx. Then 

J'(*) = 2 / [-VwfeVMfe + V k u 2 k ]dx. 
Jd 

By our assumption on V, 

J'(t) < 2bJ(t) 
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which implies 

/ u1(x,t)dx < / ul(x)dx e 2bt . 

J D J D 

Therefore if uq G L 2 (D), we conclude that Uk{x,t) increases to a finite positive limit 
u(x, t) as k — > oo, for all t and for Moreover u(-,t) G L 2 (R n ). 

Write Wk = logWfc. From ()3.1|) . one deduces 

Aw fc + |Vw fc | 2 + y fc - (w k ) t = 0. 

Therefore 

w k (x,t) = / G o (x 7 t;y,0)w k (x,0)dy 
Jd 

+ G (x,t;y,s)\Vw k {x,s)\ 2 dyds + / / G {x,t;y, s)V k (y, s)dyds. 

Jo Jd Jo Jd 

Therefore 

/ G (x,t;y,s)V k (y,s)dyds <w k (x,t) - G (x,t;y,0)w k (x,0)dy. 
'o Jd Jd 

By the monotone convergence theorem 

/ G (x,t;y, s)V(y, s)dyds < w(x,t) - / G o (x,t;y,0)w(y,0)dy 
'a Jd Jd 

< log(l + u(x,t)) - [ G (x,t;y,0)\ogu (y)dy. 

Now we take 

u (x) 



D 



1 + \x\ n 

Then, since uq G L 2 (R n ), we have 

u(-,t) G L 2 (R n ). 

By Jensen's inequality, 

log(l + u(-, t)) G L p (R n ), Vp>l. 

It is also clear that 

/ G (;t;y,0)\og(l + \y\ n )dyeLl c (R n ), Vp > 1. 
The result follows. □ 



4. APPLICATIONS TO THE NAVIER-STOKES EQUATION 

In this section, we establish a new a priori estimate for a certain quantity involving the 
velocity and vorticity of the 3 dimensional Navier-Stokes equation. 



(4.1) 



u t — Au(x, t) + u ■ Vu(x, t) + Vp = 0, 

V-M = 0, 

u(x, 0) = Mo(x) 
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2.3 



for (x,t) G IR 3 x (0, oo), where A is the standard Laplacian, a vector field u represents the 
velocity of the fluid, and a scalar field p the pressure. (The viscosity is normalized, v — 1.) 

There has been an extensive and rapidly growing literature on the equation which is 
impossible to quote extensively here. Let us just mention that weak solutions are known 
to exist due to the seminal work of Leray [T] . However it is not known if the weak solution 
is smooth everywhere. Several sufficient conditions implying smoothness of weak solutions 
have been made. See for example |F] and fSJ. In these two papers, it was shown that if the 
velocity u is in L p,q class with - + - < 1, then u is actually smooth. For more sufficiency 



results in various other spaces we refer the reader to the more recent survey paper [Ca . 
However it is only known that u G L 10 A 1( V 3 . Therefore there is a gap in between the a 
priori estimate and the sufficiency condition. 

What we will prove here is a different sufficiency condition and a priori estimate using the 
heat bounded and almost heat bounded potentials defined in the previous section. There is 
still a gap between the two conditions. However the gap seems logarithmic. More precisely, 
we have 

Theorem 4.1. Let u be a Leray-Hopf solution of the Navier-Stokes equation, which is 
classical in R 3 x (0, T). Let w be the vorticity V x u. Define the quantity 



curl{u x w) ■ w + 2|Vv/|w 


2 + l| 


2 _ 


Vw\ 


2 




w 


| 2 + 1 



Q = Q(x,t) = * '- , ' v >-(x,t) 



Then the following statements hold for any 5 G (0, T). 

(1) . The quantity Q is almost heat bounded in (R 3 x (S,T]) n {\w\ > 1}. 

(2) . u is a classical solution of the Navier-Stokes equation in R 3 x (5, T] if and only if 
Q is heat bounded in (R 3 x (8, T]) D {|u>| > 1}. 

Remark 4-1- 

The quantity Q is well defined since we assume that u is smooth for t G (0, T). The first 
term in Q is essentially the vortex stretching factor which is the hardest to control. The 
point of the theorem is that if there is blow up at time T, then the blow up just happens 
barely. 

Proof of Theorem 4-1- 

We will just prove (1) since (2) is self-evident afterward. 
We divide the proof into three steps. 

Step 1. rewriting the vortex equation in the log form. 

Let w = w(x,t) be the vortex. It is well known that \w\ 2 satisfies the following scalar 
heat equation with lower order terms 

(4.2) AH 2 - u ■ VM 2 + 2a\w\ 2 - 2\Vw\ 2 - {\w\ 2 ) t = 0. 



Here a is the vortex stretching potential given by (c.f. [Co]) 
(4.3) 



3 f dy 

a(x,t) = —P.V. / D[y,uj(x + y),uj(x)) \u(x + y,t)\ — 



wVu ■ w 



\w\ 2 
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A straightforward computation from (4.2) shows 
Am(H 2 + l)-w-Vln(M 2 + l) + 2 
Write f = \ ln(H 2 + 1). We deduce 



a 


w 


2 


\w 


2 + l 





2 , \(V\w 


2) |2 


\w 


2 + l ' (\W 


2 +l) 2 



^(ln(| W | 2 + l)) = 0. 



(4.4) 



A/ - u ■ V/ 







2 




2 + l 



T + 2|V/| 2 - 



\Vw\ 



\w\ 



- ft = o. 



Step 2. a representation formula. 

By our assumption, for t G (0,T), w and u> are classical functions and / vanishes near 
infinity. This shows, 
(4.5) " 

f(x,t) = J G (x,t;y,s)f (y)dy 



+ 



G (x,t;y,s) 



a 


w 


2 


\w 


2 + l 



«-V/ + 2|V/| : 





2 


\w 


2 + l 



(y, s)dyds. 



Step 3. Apply Jensen's inequality. 
For convenience, we write 



(4.6) 

It is clear that 



Q 



a 


w 


2 


\w 


2 + l 



M -V/ + 2|V/p 





2 


\w 


2 + l 



^-I, j9 ,MH 2 + i) + 2|v/| 

_l_ 2|V/| 2 





2 




2 + l 



w 2 + 1 





2 




2 + 1 



wVm • w — nVw • w 


+ 2IVV 




2 + 1 


2 _ 


Vw| 2 






2 


+ 1 





Following the well known vector identity, we have 
(4.7) Q = 



curliu x w) ■ w + 2 Vy 7 


ID 


2 + 1 


2 _ 




2 






2 + 1 



It is well known that w G L 2 (R 3 x R + ). Using Jensen's inequality, it is easy to show 
that /(•, t) G L p for any p > 1, in the region where > 1. Hence the quantity Q is almost 
heat bounded. □ 
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